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We show numerically that the Abelian Higgs field equations in the background 
of a four-dimensional rotating charged black string have vortex solutions. These 
solutions which have axial symmetry show that the rotating black string can support 
the Abelian Higgs field as hair. We find that one encounters with an electric field 
coupled to the Higgs scalar field for the case of rotating black string. This electric 



> 

\^ . field is due to an electric charge per unit length, which increases as the rotation 

parameter becomes larger. We also find that the vortex thickness decreases as the 
rotation parameter grows up. Finally we consider the self-gravity of the Abelian 
Higgs field and show that the effect of the vortex is to induce a deficit angle in the 



metric under consideration which decreases as the rotation parameter increases. 



I. INTRODUCTION 

The classical no hair conjecture, first proposed by Ruffini and Wheeler states that the 
only allowed characteristics of stationary black holes are those associated with the Gauss 
law, such as mass, angular momentum and electromagnetic charge jlj. Thus, for example, 
it follows that black holes cannot support external scalar fields, since there is no Gauss 
law for scalars. This fact is verified in for a scalar field minimally coupled to gravity 
and in j^] for non-minimally coupled scalar field in asymptotic flat spacetimes. However, 
recently, it is shown that there exist scalar hairy black hole if one uses a scalar potential 
with regions violating the weak energy condition j^. In recent year, we are faced with 
new black hole solutions in special gauge theories that Einstein's equation is coupled with 
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some gauge fields, and therefore this conjecture needs more investigation. For example the 
Einstein- Yang-Mills black hole possess a short ranged external non Abelian gauge field and 
is not uniquely specified by their mass, angular momentum and electromagnetic charge [51. 
Al t ho„ gh this section is unable Q, it opened np new possibility that the bl ae k holes 
may have various matter hairs. The global monopole field is an example of stable hair 
for asymptotically flat black hole Q|. Also explicit calculations have been carried out which 
verify the existence of a long range Nielson-Olesen vortex solution as a stable hair for a static 
black hole solution of Einstein's equation in four dimensions 8j. For the extreme black hole 
of Einstein- Maxwell gravity, it has been shown that flux expulsion occurs for thick strings 
thick with respect to the radius of horizon), while flux penetration occurs for thin strings 
oHlf). 11. 12^. Of course one may note that these situations fall out side of the scope of the 
classical no hair theorem due to the non trivial topology of the string configuration. 

It has recently been shown that these ideas can be extended to the case of anti de sitter 
(AdS) and de Sitter (dS) spacetimes. For asymptotically AdS spacetimes, it has been shown 
that conformally coupled scalar field can be painted as hair [13J. Another asymptotically 
AdS hairy black hole solution has been investigated in Ref. 14| . Also it was shown that there 
exist a solution to the SU(2) Einstein- Yang-Mills equations which describes a stable Yang- 
Mil, hairy blaek ho.e, that is a 8y ,nptot 1C a. ly AdS Ql Although the ,dea of Nie.son-O.esen 



vortices has been first introduced in flat spacetimes [l6( , but recently it has been extended 



to (A)dS spacetimes |17I. Il8j|. The existence of long range Nielson-Olesen vortex as hair for 
asymptotically AdS black holes has been investigated in Refs. 3, l^j] for Schwarzschild-AdS 
black hole and charged black string. 

The presence of a rotation parameter introduces significant additional difficulties for 
explicit calculations intended to investigate the existence of a long range Nielson-Olesen 
vortex solution as a hair for a stationary black hole solution j^]. In this paper we study the 
Abelian Higgs hair in a four dimensional rotating charged black string that is a stationary 
model for Einstein-Maxwell equations in cylindrical symmetry. Various features of this kind 
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23]. 



Here, we want 



of solutions of Einstein-Maxwell equations have been considered 
to investigate the influence of rotation on the vortex solution of Einstein-Maxwell-Higgs 
equations. Since an analytic solution to the Abelian Higgs field equations appears to be 
intractable, we confirm by numerical calculations that the rotating charged black string can 
be dressed by Abelian Higgs field as hair. 
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In Sec. |nj we consider the Abelian Higgs field equations in the background of a rotating 
charged black string. Section IIHI is devoted to the numerical solutions of the field equations 
for different value of the rotation parameter and winding number. In Sec. IIV| by studying 
the behavior of the Abelian Higgs field energy-momentum tensor, we find the effect of the 
vortex self-gravity on the rotating charged black string background. We append some closing 
remarks in Sec. 

II. FIELD EQUATIONS 

The Lagrangian of Einstein gravity in the presence of electromagnetic and Abelian Higgs 
fields is Q 

C = K - 2A - F^F*™ + C H , (1) 

where A = — 3// 2 is the negative cosmological constant and F^ v = d^A v — d u A^ is the field 
strength tensor associated with the electromagnetic field A^. The first three terms in Eq. 
(JIJ are the Einstein-Maxwell volume terms with negative cosmological constant, and Ch is 
the Lagrangian for Abelian Higgs system minimally coupled to gravity given as 

C H = -\{V^V^ - -±-T^ v - f - rff. (2) 

Z 107T 

The matter content of the Abelian Higgs system consists of the complex Higgs field, $, and 
a U(l) gauge field £> M with strength = d fl B u — d u B fl . Both the Higgs scalar and the 
gauge field become massive in the broken symmetry phase. They are coupled through the 
gauge covariant derivatives = V^ + ieB^, where V M is the spacetime covariant derivative. 
We use units in which 8irG = c = 1. 

The fields in Ch will be treated as "test field", i.e., their energy momentum tensor is 
supposed to yield a negligible contribution to the source of the gravitational field. Notice 
that we have two different gauge fields, F and J 7 , and each is treated in a different manner. 
It is only T that couples to the Higgs scalar field and is therefore subject to spontaneous 
symmetry breaking. The other gauge field, F, can be thought of as the free massless Maxwell 
field which apart from modifying the background geometry, its dynamic will be of little 
concern to us here. It is conventional to express the field equations in terms of real fields 
X(x^), u(x^), and P^,{x v ) defined as 

= rjXix^e™^, (3) 
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fl M 0O = i^M-vX^)]. ( 4 ) 

where X and P M are the scalar Higgs field and massive vector boson, and uj is a gauge degree 
of freedom. In terms of these new variables, the Lagrangian, Ch arid the field equations can 
be written as 

C H (X, P M ) = -|(V,XVI + X 2 P,P^) - -^F'^ - ^{X 2 - 1), (5) 
V^VI - XP/" - 4£r/ 2 X(X 2 - 1) = 0, (6) 

VfiF '^ _ ^e 2 rfP v X 2 = 0, (7) 
3 

G^u — pQiJiv = + 99 , (8) 

where F ,fiu = V^P 1 * — V^P M is the field strength of the corresponding gauge field P M , and 
and F^ v % " s are the stress energy tensors of the electromagnetic and Higgs fields given 

by 

^jtv = 2P Al P 7/ — —F g^ u , (9) 

T Higgs = rf Vii XV v X + 7 1 2 X 2 P I1 P V + -^FvrFZ + ^„£tf. (10) 

Note that the real field uj is not itself a physical quantity. Superficially it appears not to 
contain any physical information. However, if uj is not single valued this is no longer the 
case, and the resultant solutions are referred to as vortex solutions jlfil ]. In the presence of 
vortex, the requirement that the $ field be single valued implies that the line integral of uj 
over any closed loop is ±27rN, where N is the winding number of the vortex. In this case 
the flux of the field T passing through such a closed loop is quantized with quanta 27r/e. 

It must be noted that no exact solution are known for the coupled Einstein-Higgs equa- 
tions ©-(0) even in situations having both no electromagnetic charge and no horizon. Thus 
we should try to solve these coupled differential equations approximately. In the first order 
approximation, we solve Eqs. (jHJ) and (J2J) in the background of charged rotating black string 
and then we will carry out numerical calculations to solve Eq. (jHJ) . The metric of the charged 



rotating black string can be written as 
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24] 



ds 2 = -r (Edt - ad(f>) 2 + — (adt - El 2 d(f) + — + — dz 2 , (11) 

I Li 

where 

,2 U7 \2j2 



r 2 bl X 2 l 2 . \l 



[2 r r 2 



A, = -- (aSj - aS*) , E=^lTW 2 - (12) 
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b and A are the constant parameters of the metric which are related to the mass and charge 
per unit length of the black string by M = 6/4 and Q — A/2, and a is the rotation parameter. 
It is remarkable to note that the metric (fTTj) reduces to AdS spacetime as t goes to infinity 
in contrast to the Kerr-AdS metric, for which the surface at infinity is also rotating. The 
horizon of this black string has cylindrical symmetry that has two inner and outer horizon 
at r_ and r+, provided that b is greater than b crit given by 



J crit 



4 x 3- 3 / 4 A 3 / 2 . (13) 



In the case that b = b cr i t , we will have an extreme black string. 

We seek a cylindrically symmetric solution for the Higgs field equations (jSJ) and ((7j) in 
the background of a charged rotating black string. Thus, we assume that the fields X and 
Pfj, are functions of r. For the case of vanishing rotation parameter, one can choose the 
gauge field as P M (r) = (0, 0, Np(r), 0). The field equations (jSJ and (JZj) reduce in that case to 
two equations for the two unknown functions X(r) and P(r). In the present case the field 
equations (J5J) and (JJJ) reduce to three equations and one may therefore use the gauge choice 

P,(r) = (S(r)ANP(r),0). (14) 

The field equations fl^J) and (JJJ) reduce to 

r 2 r X" + (4r 2 r 3 - bl) X' - 4 r 2 X (X 2 - l) - X (aS + ENP) 2 - 

-r 2 r -1 X (ES + Nal~ 2 P) 2 = 0, (15) 

N{rH 2 TE 2 - r b a 2 r 2 )P" + iV(2r 4 + 6/ 3 rS 2 - 2 A 2 / 4 S 2 )P' + (Z 2 H 2 - rVr 2 !^ 1 ) x 

Nar 3 X 2 P + l 2 aE(X 2 l 2 - Mr) (rS" - S' + a rX^S) - aEX 2 l 4 S' = 0, (16) 

r 3 (~V -Ta 2 ) 5" + [2r 4 - a 2 (Wr -2X 2 l 2 )]S' + aX 2 r 3 T- 1 (~V - la 2 ) S 

—NE a[{rP" - P' + ar X^P) (A 2 / 2 - Mr) - \ 2 l 2 P'} = 0, (17) 

where a = 47re 2 /£ and the prime denotes a derivative with respect to r. It is worthwhile to 
mention that even in the pure flat or (anti) de Sitter spacetimes no exact analytic solutions 
are known for equations (jHJ) and Q. For asymptotically AdS spacetimes, one of us have 
showed that the Abelian Higgs equations of motion in the background of charged black string 
spacetime have vortex solution 2^|. Here we want to investigate the influence of rotation 
parameter on the vortex solutions. 
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FIG. 1: Effiggg x 10 3 versus r for N = 1, a = (touch the horizontal axis), 0.25 (dotted), 0.5 
(solid) and 0.7 (bold). 

III. NUMERICAL SOLUTIONS 



The main difference between this rotating black string and the non rotating case consid- 
ered in Ref. j^J is the following. In the static case we had only two unknown functions 
X(r) and P(r), but here we encounter with three unknown functions X(r), P(r) and S(r). 
The latter function S(r) which is the time component of creates a radial electric-type 
field, F tr that goes to zero as the rotation parameter goes to zero. It should be noted that 
even in the case of b = A = 0, no exact analytic solutions are known for Eqs. (|15 p -()17j ) . 
So, here we seek the existence of vortex solutions for these coupled differential equations 
numerically. First, we must take appropriate boundary conditions. Since at a large distance 
from the horizon the metric (jllj) reduces to AdS spacetime, we demand that our solutions 
go to the solutions of the vortex equations in AdS spacetime given in fl^ |. This requires 
that we demand (X — > 1, P — > 0) as r goes to infinity and (X = 0, P = 1) on the horizon. 
For consistency with the non-rotating case j^, we take S = on the horizon. Also one 
may note that the electric field, F' tr which is proportional to 5" should be zero as r goes to 
infinity. We employ a grid of points r^, where r« goes from r# to some large value of r (r^) 
which is much greater than r#. We rewrite Eqs. fll5|) - (fT7j) in a finite deference language 
and use the over relaxation method to calculate the numerical solutions of X(r), P(r) 
and S(r) for different values of the rotation parameter and winding number. The numerical 
results of calculations are shown in Figs EE 

First, we investigate the influence of rotation parameter on the solutions of field equations 
(|15 )l -p7 |l . We carry out all the calculation for I — 1, A = 0.2 and b = 0.3 for which the radius 



FIG. 2: H Higgs versus r for N = 1, a = (solid) and 0.7 (bold). 
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FIG. 3: E Hi g gs x 10 3 versus r for JV = 1 (dotted), 3 (solid), and 5 (bold). 

of horizon is r H = 0.6173. Figures [T] and El show the behavior of the electric, E Higgs = F t ' r , 
and the magnetic, Hmggs = FL,, fields associated with the field P M respectively for different 
values of the rotation parameters. We use the subscript "Higgs" for these electromagnetic 
field to emphasis that they are coupled with the Higgs scalar field $. As one can see in 
Fig. HI Ejiiggs is zero for a = 0, and becomes larger as a increases. The magnetic field 
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FIG. 4: H Higgs versus r for = 1 (dotted), 3 (solid), and 5 (bold). 
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FIG. 5: X(r) versus r for N = 1 (dotted), 3 (solid), and 5 (bold). 



Huiggs is plotted in Fig. El for different values of angular momentum. This figure shows that 
the variation of HHi ggs with respect to the rotation parameter is very slow. Overall, these 
figures show that the vortex thickness decreases as the rotation parameter increases. As we 
mentioned, in the case of non vanishing rotation parameter we encounter with the electric 
field Euiggs- One may compute the source of this electric field through the use of Gauss 
law numerically. It is notable that the computation of this electric type charge for Em g gs 
shows that it increases as a becomes larger. This is analogous to the rotating solutions of 
Einstein-Maxwell equation discussed in the context of cosmic string theory for which the 
electric charge of the string is proportional to the rotation parameter of the string (26 1. 

Next we investigate the influence of the winding number N on the solutions of field 
equations (fTH|) - (fT7j) for the case of rotating charged black string. The results for a = 0.5 are 
shown in Figs. 01 El and El for different values_of N. As in the case of asymptotically flat, 



and 



17 



la . ll9L |20| , increasing the winding 



dS, and AdS spacetimes considered in Refs. 
number yields a greater vortex thickness. 

The effects of charge or mass parameters on the solutions of field equations (|15|) -(|17 J) 
for the case of rotating charged black string are the same as the nonrotating case, which 
is considered in j^j. The extreme black string needs more investigation. Bonjour et.al. 



12j have shown that vortices of size smaller than a certain fraction of the radius of the 
Reissner-Nordstrom black hole will definitely pierce the horizon, whereas vortices thicker 
than a certain lower bound will instead wrap the black hole. We will consider the extreme 



case in future. 
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IV. ASYMPTOTIC BEHAVIOR OF THE SOLUTIONS OF 
EINSTEIN-MAXWELL-HIGGS EQUATION 

In previous section we found the solutions of Higgs field equation in the background of 
charged rotating black string. Here we want to solve the coupled Einstein-Maxwell-Higgs 
differential equation This is a formidable problem even for flat or AdS spacetimes, 

and no exact solutions have been found for these spacetimes yet. Indeed, besides the elec- 
tromagnetic stress energy tensor, the energy-momentum tensor of the Higgs field is also a 
source for Einstein equation (jHJ. However, some physical results can be obtained by making 
some approximations. First, we assume that the thickness of the skin covering the black 
string is much smaller than all the other relevant length scales. Second, we assume that 
the gravitational effects of the Higgs field are weak enough so that the linearized Einstein- 
Abelian Higgs differential equations are applicable. We choose ~ g$ +sg$, where g$ 
is the rotating charged black string metric in the absence of the Higgs field and g$ is the 
first order correction to the metric. Employing the two assumptions concerning the thick- 
ness of the vortex and its weak gravitational field, the first order approximation to Einstein 
equation (jHJ) can be written as: 

G$-p$ = T$\ (18) 

where G^} is the first order correction to the Einstein tensor due to g$ and iffl is the 
energy-momentum tensor of the Higgs field in the rotating charged black string background 
metric with components: 

T t m (r) = {-H 4 (r 2 E 2 -Ta 2 )S' 2 - ar 2 fV 2 X fl - N 2 T (eYT - a 2 r 2 ) P' 2 

-a N 2 X 2 fYP 2 - aa 2 (T I 2 - r 2 ) (N 2 P 2 - S 2 l 2 ) X 2 - arY X 2 S 2 

-2ar 2 l 4 T (X 2 - l) 2 }/2r 2 / 4 r, 
T/ (0) (r) = { H 4 (-r a 2 + r 2 E 2 ) S' 2 - ar 2 l 4 T 2 X' 2 + N 2 T (T E 2 l 4 - a 2 r 2 ) P' 2 

+ al A X 2 (N 2 TP 2 + r 2 S 2 ) + a X 2 a 2 (T I 2 - r 2 ) (N 2 P 2 - l 2 S 2 ) 

-2aTr 2 l A (X 2 -l) 2 }/2r¥r, 
r; (0) (r) = {2N 2 l 2 TaE (T I 2 - r 2 ) P' 2 + 2H 4 iV (r a 2 - r 2 S 2 ) P'S' - r 2 l 2 PS 

+2l 2 NX 2 a (Pa (T I 2 - r 2 ) (aS + E NP)) } /2r 2 / 4 r, 

r;(°)(r) = {n 4 (r a 2 - hV) s' 2 + arYr 2 x' 2 + n 2 t (teY - a\ 2 ) p' 2 
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+2a2/ 2 AT (H 2 - r 2 ) P'S' - 2ar 2 l 4 T (X 2 - l) 2 - a X 2 a 2 x 
(r I 2 - r 2 ) (N 2 P 2 + l 2 S 2 ) - 2a l 2 X 2 Na~ SP (T I 2 - r 2 ) 
- al 4 X 2 (TN 2 P 2 - r 2 S 2 ) } /2r 2 / 4 I\ 
T<°\r) = {Tl 4 (-Ta 2 + r 2 E 2 )S' 2 - ar 2 l 4 T 2 X' 2 - N 2 T(Tl 4 E 2 -a 2 r 2 )P' 2 
-2al 2 ~NT (VI 2 - r 2 ) P'S' - 2ar 2 l 4 T (X 2 - l) 2 - alV x 
(r I 2 - r 2 ) (N 2 P 2 + l 2 S 2 ) - 2a/ 2 X 2 Na~ SP (T I 2 - r 2 ) 
-a l 4 X 2 (N 2 TP 2 - r 2 S 2 ) } /2r 2 l 4 T, (19) 

where X, P, and S are the solutions of the Abelian Higgs system. As we mentioned in the 
last section, the vortex thickness decreases as the rotation parameter increases. This fact is 
more clear in Fig. |H1 which shows for various values of a. 

For convenience, we use the following form, which has cylindrical symmetry, of the metric 

ds 2 = -A(r)dt 2 + B(r)dr 2 + C(r)dtd(f) + D(r)d(j) 2 + E(r)dz 2 . (20) 

In order to solve numerically Eq. (JIEJ), it is better to write the metric function A(r) 
to E{r) as 





A(r) 


= A (r)[l + eA(r)], 




B(r) 


= B (r)[l + eB(r)], 




C{r) 


= C (r)[l + eC(r)}, 


D(r) 


= D (r)[l+eD(r)), 




E(r) 


= E (r)[l + eE(r)\, 
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FIG. 6: T t versus r for a = (solid) and a = 0.7 (bold). 
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FIG. 7: T/ (0) = T/ (0) (broken), T/ (0) (solid), T * {O) (bold) and T; (0) (thick-bold), versus r for 
iV = 1, a = 0.5;. 
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FIG. 8: A(r), B{r) and S(r) touch the horizontal axis, C(r) (solid) and D[r) (bold). 

where A Q (r) = TE 2 - r 2 a 2 l~\ B {r) = F' 1 , C = 2aE(T - r 2 r 2 ), D (r) = r 2 E 2 - Fa 2 , 
Eq(v) = r 2 l~ 2 , yielding the metric of the stationary rotating charged black string in four 
dimensions. The Einstein equations (fT%J) in terms of the functions A(r) to E(r) are given 
in the appendix. Here we want to obtain the behavior of these functions for large values 
of the coordinate r. As one can see from Fig. [7[ the components of the energy-momentum 
tensor rapidly go to zero outside the skin, so the situation is the same as what happened in 
the static black string spacetime considered in (^(J. One can solve the linearized Einstein 
equation for large values of r numerically. The results which are displayed in Fig. |H1 show 
that A(r) = B(r) = E(r) = 0, and 2C(r) = D[r) = 2. Hence the metric ()2()|l can be written 
as 

ds 2 = -A (r)dt 2 + B (r)dr 2 + (1 + e)C (r)dtd<p + (1 + 2e)D (r)d(p 2 + E (r)dz 2 . (22) 

It is worthwhile to mention that the metric ()22j) is the first order solution in e of the Einstein- 
Maxwell-Higgs equations far from the thin string. Of course, one may note that the metric 
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(122)1 is the first order approximation of the following metric 

ds 2 = -r (Edt - a(3d<j)) 2 + — (adt - El 2 pd<p) + — + —dz 2 , (23) 

which is the exact solution of Einstein-Maxwell gravity. In Eq. (|23)1. (3 is defined as (3 = 1+e. 
The above metric describes a stationary rotating charged black string with a deficit angle 
A = 2ns. The size of deficit angle A is proportional to 2n J rT^dr Numerical com- 
putation shows that the absolute values of this integral decreases as the rotation parameter 
increases, which can also been from Fig. |H1 So, using a physical Lagrangian based model, we 
have established that the presence of the Higgs field induces a deficit angle in the rotating 
charged black string metric which decreases as the rotation parameter increases. 



V. CONCLUSION 



The effect of a vortex on pure AdS spacetime is to create a deficit angle in the metric 
in the thin vortex approximation. This fact has been extended to the charged rotating 
spacetimes. This is done by including the self-gravity of the thin vortex in the rotating 
charged black string background in the first order approximation discussed in Sec. II VI As 
in the case of pure AdS |l2(, Schwarzschild-AdS jlSj, Kerr-AdS, and Reissner-Nordstrom- 
AdS [21] spacetimes, we found out that the effect of a thin vortex on the stationary charged 
black string is to create a deficit angle in the metric which decreases as the rotation parameter 
increases. 

We considered the Abelian Higgs field in the background of a stationary rotating charged 
black string. We obtained the numerical solutions for various values of rotation parameter 
and found that for a fixed horizon radius, by increasing the rotation parameter the vortex 
thickness decrease very slowly. Also the numerical solutions for various values of winding 
number were obtained. These solutions shows that the vortex thickness increases as winding 
number increases. These situations hold also in the case of the extremal black string. 

The main difference between the case of Abelian Higgs field in the background of static 
black string considered in , and this paper is that the time component of the gauge field 
coupled to the Higgs scalar field is not zero for non zero rotation parameter. Indeed, we 
found that for the case of rotating black string, there exist an electric field coupled to the 
Higgs scalar field. This electric field increases as the rotation parameter becomes larger. 
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Numerical calculations show that the electric charge which creates this electric field grows 
up as the rotation parameter increases. This is analogous to the results that Dias and Lemos 
have found recently for the magnetic rotating string [26]. They showed that the charge per 
unit length of a rotating string in the Einstein-Maxwell gravity increases as the rotation 
parameter becomes larger, and we found that the electric charge of the field F' coupled to 
the Higgs field has the same feature. 

Other related problems such as a study of extreme black string in the presence of Abelian 
Higgs field, a study of the non-Abelian vortex solution in asymptotically AdS spacetimes, 
or other kinds of fields coupled with gravity remain to be carried out. Work on these 
problems is in progress. 



[1] R. Ruffini and J. A. Wheeler. Phys. Today 24, 30 (1971). 
[2] D. Sudarsky, Class. Quantum Grav. 12, 579 (1995). 

[3] A. Saa, J. Math. Phys. 37, 2346 (1996); A. Saa, Phys. Rev. D 53, 7377 (1996); A. E. Mayo 

and J. D. Bekenstein, Phys. Rev. D 54, 5059 (1996). 
[4] U. Nucamendi and M. Salgado, Phys. Rev. D. 68, 04026 (2003). 

[5] M. S. Volkov and D. V. Gal'tsov, Y. Fiz. 51, 1171 (1990) [Sov. J. Nucl. Phys. 51, 747 (1990)]; 
P. Bizon, Phys. Rev. Lett. 64, 2844 (1990); H. P. Kunzle and A. K. Masoud-ul-Alam, J. Math. 
Phys. 31, 928 (1990). 

[6] N. Straumann and Z. H. Zhou, phys. Lett. B 243,33 (1990); Z. H. Zhou and N. Straumann, 
Nucl. Phys. B 360, 180 (1991); P. Bizon, Phys. Lett. B 259, 53 (1991); M. S. Volkov, et al, 
Phys. Lett. B 349, 438 (1995); O. Brodbeck, N. Straumann, J. Math. Phys. 37, 1414 (1996). 

[7] H. Watab and T. Torii, JCAP 0402, 001 (2004). 

[8] A. Achucarro, R. Gregory, and K. Kuijken, Phys. Rev. D 52, 5729 (1995). 
[9] A. Chamblin, J. M. A. Ashbourn-Chamblin, R. Emparan, and A. Sornborger, Phys. Rev. D 
58, 124014 (1998). 

[10] A. Chamblin, J. M. A. Ashbourn-Chamblin, R. Emparan, and A. Sornborger, Phys. Rev. Lett. 
80, 4378 (1998). 

[11] F. Bonjour and R. Gregory, Phys. Rev. Lett. 81, 5034 (1998). 



14 



[12] F. Bonjour, R. Emparan, and R. Gregory, Phys. Rev. D. 59, 84022 (1999). 

[13] E. Winstanley, Found. Phys. 33, 111 (2002). 

[14] T. Torii, K. Maeda, and M. Narita, Phys. Rev. D 64, 042116 (2001). 

[15] E. Winstanley, Class. Quantum Grav. 16, 1963 (1999). 

[16] H. B. Nielsen and P. Olesen, Nucl. Phys. B 61, 45 (1973). 

[17] M. H. Dehghani, A. M. Ghezelbash, and R. B. Mann, Nucl. Phys. B 625, 389 (2002). 

[18] A. M. Ghezelbash and R. B. Mann, Phys. Lett B 537, 329 (2002). 

[19] M. H. Dehghani, A. M. Ghezelbash, and R. B. Mann, Phys. Rev. D 65, 044010 (2002). 

[20] M. H. Dehghani and T. Jalali, Phys. Rev. D 66, 124014 (2002). 

[21] A. M. Ghezelbash and R. B. Mann, Phys. Rev. D 65, 124022 (2002). 

[22] M. H. Dehghani, Phys. Rev. D 66, 044006 (2002). 

[23] M. H. Dehghani, A. Khodam-Mohammadi, Phys. Rev. D 67, 084006 (2003). 

[24] J. P. S. Lemos and V. T. Zanchin, Phys. Rev. D 54, 3840 (1996). 

[25] W. H. Press, S. A. Teukolsky, W. T. Vetterling and B. P. Flannery, Numerical Recipes in 

FORTRAN (Cambridge University Press, Cambridge, England, 1992). 

[26] O. J. C. Dias and J. P. S. Lemos, Class. Quantum Grav. 19, 2265 (2002). 

APPENDIX 

The Einstein equation mentioned in Sec. (jlVjl in terms of the metric functions A(r) to 
E(r) are 

E 2 Fi(A + D-2C)- 4A Wr 2 (r S 2 / 4 - r 2 a 2 )D + 4A 2 / 8 ~ 2 r 2 (r a 2 - r 2 ~ 2 ) A - rH 2 T 2 x 
{4H 4 S 2 r + (bl 3 - 4r 3 a 2 ) } B' - E 2 a 2 l s rT 2 {2 b 2 r 2 + A 2 (4A 2 / 2 - 6 bit)} A' - 12Z 4 r 6 T 2 B 
- Z 2 r 2 l 2 TF2 (A' + D'-2C') + / 2 r 2 r 4 {4r/ 4 r S 2 + (4r 3 - bl 3 ) (a 2 - I 2 ) } (E' + D') 
+2TW {E" + D") + T 2 l A a 2 {2-E?l A [b 2 r + \ 2 {2\ 2 l 2 r- 1 - 36/)] + r 3 (4A 2 / 2 - 3blr)}D' 
-2~ 2 a 2 \ 2 l 6 rT 2 (Trl 2 - r 3 - bl 3 ) (D" - C") = 4r 6 / 6 r 2 ^ (0) , 

8r 2 ( Yl 2 - r 2 )A 2 H fa 2 D + 8 A 2 S Z 6 r 2 a{Z 2 r 2 - (Tl 2 - r 2 )a 2 }C + aH 3 F 3 (A + D - 2 C) + 

r 4 / 4 aSr 2 (4H 2 r - 4r 3 + bl 3 ) (B' — E' — 2D') - r 2 l 2 aZTF 5 (A' - D') - 2r 2 / 2 Harr 4 x 

(C" _ D ') + 2 / 4 r 2 S a {/V 2 (S 2 - 1) r (A 2 / - for) + S 2 r 6 (r 2 - Tl 2 )} (D" - C") = Ar G l 6 T 2 T^ \ 

F 6 (A + D-2C) + 8 a 2 \ 2 l 6 T 2 {Tl 2 E 2 - r 2 (-1 + S 2 ) } (C - A) + Ar 2 l 6 \ 2 T 2 ( E 2 l 2 + a 2 ) A+ 
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r 2 rn 4 H 2 F 7 (A' + D' - 2 C) - AE 2 a 2 l 6 rT 2 {b 2 rl 2 + A 2 (2 Trl 2 - bl 3 - 2r 3 )}C + l A r 2 T 2 F 8 A' 
+r 4 / 4 r 2 (4r 3 - bl 3 )E' + r 4 / 4 r 2 {4r a 2 T - H 2 (4r 3 - bl 3 )} {B' - E') + 2r 6 / 6 T 3 (E" + A") 
-2E 2 a 2 X 2 l 6 rT 2 (Tl 2 r - r 3 - bl 3 ) (A" - C") - 12l A r 6 T 2 B = 4r 6 l G T 2 T^°\ 

E 2 F 9 (A + D - 2 C) - 8 r 2 A 2 S 2 / 6 rC + 4r 2 A 2 / 4 r ( E 2 l 2 - a 2 ) D + r 4 / 2 r (4r 3 - bl 3 ) E' + r 4 l 2 T 
x{4r a 2 T - ~ 2 (4r 3 - bl 3 )}D' + r 4 r{4H 4 S 2 r - a 2 (4r 3 - bl 3 )}A' + S 2 a 2 6/ 3 r 2 r(H 2 - r 2 ) 
x (A' + D' — 2C) - l2rH 2 YB = 4r 6 / 4 rT r r(0) , 

S 2 rio (A + D-2C)- 16r 2 A 2 / 6 r 2 ( ~ 2 l 2 - a 2 ) D + (32 A 2 ~ 2 /Vr 2 )C - 4/ 4 r 4 r 2 (4r 3 - bl 3 )B' + 
4/ 4 r 2 r 2 {8r 3 / 2 V + (4A 2 / 2 - 3blr) (a 2 - 2l 2 ) }D' + 4/Vr 2 {8 r 4 + (3 E 2 - 2) bl 2 r - 4A 2 / 4 S 2 }A' 
+4/ 2 r 2 r H 2 Fn (A' + D' — 2C) - 8 H 2 a 2 A 2 / 6 rr 2 (rl 2 T - r 3 - bl 3 ) (A" + D" - 2C") + 8 r 6 / 6 r 3 
x [A" + D") - 48/ 4 r 6 r 2 5 = 16r 6 / 6 r 2 T 2 2(0) , 

where T^ ^ is the energy momentum tensor of the Higgs field given in Eqs. (|19|) and the 
functions F %'s are: 

Fx = H 4 (32 r 3 /V - r 2 / 10 6 2 + rW - 32T l 2 r 8 - 8 r 3 Z 9 r& + 8r 2 / 7 r 3 6 + 16 r 10 - 16 T 4 l 8 r 2 - 8r 7 bl 3 
+8 T l 5 r 5 b) + ~ 2 (r¥V + 12 r¥V 2 - 28 r¥r 4 - 12 r 2 / 4 r 6 + 60 r l 2 r 8 - 32 r 10 + 16 r 7 bl 3 
-12T l 5 r 5 b + 4 T 3 l 9 rb - 8 T 2 l 7 r 3 b - 2 rW) + 4 T 4 l 8 r 2 + r% 2 l 6 - 8 r 7 bl 3 + 4 T 3 l 9 rb + 16 r 10 
+4 T / 5 r 5 6 + 12 T 2 / 4 r 6 - 28 T l 2 r 8 - 4 T 3 l 6 r\ 

F 2 = 5 4 (4r 7 - 4 r¥r - 6/ 3 r 4 - 6/ 7 r 2 + 2 W 5 r r 2 + 12r 2 Z 4 r 3 - 12 F l 2 r 5 ) + E 2 {2 r A bl 3 - 8 r 7 - bl 7 T 2 
+ 16r 5 r I 2 - 8 r 3 r 2 / 4 - r 2 6/ 5 r) - r A bl 3 + 4 r 7 + 2 T 2 / 7 6 - Y l 5 r 2 b + 4 r 3 / 6 r - 4 T 2 / 4 r 3 - 4 T l 2 r 5 , 
F 3 = H 2 (32r l 2 r 8 - 32r 3 / 6 r 4 + 8 r 7 bl 3 - r 4 6 2 / 6 - 8 T 2 l 7 r 3 b - %Y fr% + 8 T 3 l 9 rb + T 2 l w b 2 - 
16 r 10 + 16 T 4 / 8 r 2 ) + 8 T l 5 r 5 b + 16 r 10 + 32 r¥r 4 + 8 T 2 l 7 r 3 b - 8 T 3 l 9 rb - 16 T 4 / 8 r 2 - 32 T l 2 r 8 
+r VZ 6 -8r 7 W 3 -6V°r 2 , 

F 4 = 5 4 (4r 7 - 6/ 7 r 2 - bl 3 r 4 + 12 T 2 / 4 r 3 - 12 T l 2 r 5 - 4 T 3 l 6 r + 2 bl 5 T r 2 ) + H 2 (8 r 5 T I 2 - bl 7 T 2 

F 5 = E\4T 3 l 6 r - 4r 7 + 6/ 7 r 2 - 2 bl 5 T r 2 - 12 T 2 / 4 r 3 + 12 T l 2 r 5 + bl 3 r 4 ) + E 2 (2r 2 bl 5 T - W 7 r 2 

+4r 7 + 12 r 3 r 2 / 4 - 12 r 5 T I 2 - 4 r¥r - r 4 6/ 3 ) - T l 5 r 2 b - 4 T 2 / 4 r 3 + 4 T l 2 r 5 , 
T 6 = S 6 (6 2 / 10 r 2 + 16r W - 32r 4 r 3 / 6 + 32 r 8 r I 2 + 8r 7 bl 3 + 8 r 3 Z 9 r& - b 2 l 6 r* - 8 r 5 T H - 16 r 10 



-8r 3 r 2 / 7 6) + H 4 rV/ 6 + 2 H 4 (46/V - 18 r 8 T I 2 + 6bl 7 T 2 r 3 - 8 bl 9 T 3 r + 30r 4 r 3 / 6 + 6 bl 5 T r 5 
-4r 6 r 2 / 4 + 8r 10 - b 2 l 10 T 2 + 16rW) + 4 S 2 (r 6 r 2 / 4 - 6r 4 r 3 / 6 + r 8 T I 2 - r 5 Tl 5 b + 2r 3 / 9 r6 
+4rW + r 2 / 10 6 2 /4) + 4( r 2 /V - r 3 / 6 r 4 - r 2 / 7 r 3 6), 
F 7 = 5 4 (11 rbl 5 X 4 - 10 r WA 2 - 4 A 6 / 6 + 3 r 3 b 3 l 3 ) + S 2 (3 r 6 b 2 - 7 r 5 bl\ 2 - 6 r 3 b 3 l 3 - 22 rbl 5 X 4 

+8 A 6 / 6 + 4 r 4 A 4 / 2 + 20 r 2 b 2 l 4 \ 2 ) - 3 r% 2 + 7 r 5 6/A 2 - 4 r 4 A 4 / 2 + 1 1 rbl 5 A 4 + 3 r 3 6 3 / 3 - 4 A 6 / 6 
-lOrWA 2 , 

T 8 = 4 T l 2 r 3 + 2 Z 3 r 2 H 2 6 + 8 r 3 S 2 r I 2 - 4 Z 3 H 4 6r 2 + 8 S 4 r 5 + 4 / 5 S 4 6r - 4 / 5 S 2 6r + 8 S 4 r 2 / 4 r - 

16 S 4 r l 2 r 3 - 8 S 2 r 2 / 4 r + 4 r 5 - r 2 bl 3 , 
F 9 = H 2 (12 rW - 24 r 6 r Z 2 - 4 l 7 bT 2 r - Z 8 6 2 r + 4 r 3 6/ 5 r + 12 r 8 ) - 16 r 4 r 2 / 4 + 28 r 6 T I 2 + b 2 l 8 T 

-8r 3 Tl 5 b + 4bl 7 T 2 r - 12 r 8 , 
Fio = 24 r 3 Z 9 r& - 16 T 2 l 7 r 3 b + H 2 (76 T 3 fr A - 40 T 4 / 8 r 2 + 20 T 2 l 7 r 3 b - 48 T 2 / 4 r 6 + 28 V l 2 r 8 - 16 r K 
-24 T 3 / 9 r6 - 2 r 2 / 10 6 2 - 4 T l 5 r 5 b - r A b 2 f + 8 r 7 bl 3 ) + 40 T 4 / 8 r 2 - 72 T 3 / 6 r 4 + 44 T 2 / 4 r 6 + 
16 r 10 + 4 T l 5 r 5 b - 28 T l 2 r 8 + r 4 6 2 / 6 - 8 r 7 bl 3 + 2 6 2 Z 10 r 2 , 
Fn = H 2 (4 r 7 - r%l 3 - 28 r 3 r 2 / 4 - 5 r 2 bl 5 T + 6 6/ 7 r 2 + 8 r 5 r Z 2 + 16 T 3 l 6 r) + r 4 6/ 3 - 8 F l 2 r 5 + 
5Tl 5 r 2 b - 4 r 7 - 6 T 2 l 7 b + 28 T 2 l 4 r 3 - 16 T 3 / 6 r. 



